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Abstract 

The field-induced transition in one-dimensional S = 1 Heisenberg antiferro- 
magnet with single-ion anisotropy in the presence of a transverse magnetic 
field is obtained on the basis of the Schwinger boson mean-field theory. The 
behaviors of the specific heat and susceptibility as functions of temperature 
as well as the applied transverse field are explored, which are found to be 
different from the results obtained under a longitudinal field. The anomalies 
of the specific heat at low temperatures, which might be an indicative of a 
field-induced transition from a Luttinger liquid phase to an ordered phase, 
are explicitly uncovered under the transverse field. A schematic phase dia- 
gram is proposed. The theoretical results are compared with experimental 
observations. 

PACS Numbers: 75.50.Ee, 75.30.Gw, 75.40.Cx 
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I. INTRODUCTION 



The study on one-dimensional (ID) magnetic systems began in 1930's. 1 Although sev- 
eral decades passed, the low-dimensional antiferromagnets have still received much renewed 
attention owing to both experimental and theoretical efforts in the past few years. On the 
theoretical aspect, most recent works associated with methods such as finite-size scaling, 3 
numerical calculations, 4,5 Monte Carlo methods, 6 analyses of an exactly solvable model, 7 
etc., have been devoted to verifying Haldane's conjecture 2 that the uniform integer spin 
chain with Heisenberg antiferromagnetic interactions is massive whereas the half-integer 
spin chain is massless, and much progress towards understanding Haldane's scenario has 
been made so far. 

On the other hand, Haldane's conjecture has been confirmed experimentally in a number 
of ID antiferromagnets with spin integer. Quite recently, Honda et al 8 ' 9 reported an anomaly 
of heat capacity on a single crystal sample of the S = 1 quasi- ID Heisenberg antiferromagnet 
(HAF) Ni(C 5 H u N 2 ) 2 N 3 (PF 6 ) (briefly NDMAP) in applied magnetic fields, and identified 
such an anomaly as an indication of a field-induced magnetic long-range ordering. They 
found an anisotropy in the susceptibility which can be explained as due to the single-ion 
anisotropy of Ni 2+ . The model Hamiltonian they adopted is expressed as 

h = j^Si-Sj + d J2(sn 2 ~ 1Mb E Si • i • h (i) 

(ij) i i 

where J is the exchange integral (for NDMAP J/ks ~ 30K), D is the single-ion anisotropy 
constant, \ib is the Bohr magneton, and g is the g tensor with elements g± and g\\ corre- 
sponding to the g values perpendicular and parallel to the chain c axis, respectively. 

Motivated by this nice experiment, we have recently discussed such an ID system in a 
longitudinal magnetic field. 10 As this magnetic system is expected to reveal quite different 
behaviors in longitudinal and transverse applied fields, as manifested by experiments 8 ' 9 , it 
would be interesting to pay attention to the case in a transverse magnetic field. It is thus 
the purpose of this present article to report thermodynamic behaviors of this system in the 
presence of a transverse magnetic field by means of the modified Schwinger boson mean-field 
theory (SBMFT). As is well known, the SBMFT works well for the ID HAF with integer 
spins, which can include the effects of quantum fluctuations self-consistently in the large- A 
saddle point, where A is the quasiparticle degeneracy. The Schwinger boson approach is 
actually a large- A formulation, and the large- A limit of Hamiltonian is taken with k — n^/N 
fixed 12 ' 13 ' 14 with the constraint b + b = at each site for Schwinger bosons. For large values 
of k the system is magnetically ordered in high dimensions, while the quantum disordered 
state appears at small k. In the presence of a magnetic field one may expect the degeneracy 
to be lifted and A could be smaller than in the absence of a magnetic field. By noting that 
the degeneracy is for Schwinger boson quasiparticles, the SBMFT still works even if the 
applied magnetic field lifts some degeneracies, because even for small A it is known that the 
SBMFT can produce fair results for the ID HAF with integer spin (see, e.g. Ref.[12]). Our 
so-obtained results would thus be reliable at the mean-field level. 

The outline of this paper is as follows. In next section, the thermodynamic quanti- 
ties such as the specific heat and the transverse susceptibility will be calculated based on 
the framework of the SBMFT. The detail formalism of SBMFT for this present system is 
collected in the Appendix. Finally, a summary and discussion will be presented. 
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II. THERMODYNAMIC PROPERTIES AT FINITE TEMPERATURES 



Based on the self-consistent equations developed in the Appendix, we may obtain the 
static uniform transverse susceptibility (along x direction) per site by 

x/2N = ^£{/? 1 5^( fc) {Mfc) + 1M*) + \n P {k) + l]n,(k)} 

+ (1 ~ ' U 2 2 ^ )] '\ n a (k) + np (k) + 1)}, (2) 

where 

U(k)=r ]lk - fiB ^ x , n a(fS) (k) = * — . (3) 

The specific heat per site, C/2N, can be easily obtained by ^p- = ^ff , where E is the 
internal energy. For convenience in numerical calculations, we shall set J = 1 as energy 
scale, g± = 2.17 (from Ref. [8]) and ks — 1 hereafter. 



A. Specific Heat 



The temperature dependence of the specific heat (C v ) for different fields at D = is 
shown as in Fig. 1. One may see that the specific heat increases with increasing field at 
low temperature, whereas it decreases with increasing field at high temperature. At a given 
magnetic field < h x < 0.9, there is an anomaly clearly observed in the curve of C v at a 
temperature T c . If the transverse field is taken off, the anomaly disappears. Such an anomaly 
might be an indicative of field-induced transition. Recall that in the presence of a longitu- 
dinal field, the system shows no any anomaly in the curve of C v observed. 10 In this sense, 
this kind of field-induced transition appears only when a transverse magnetic field is applied 
to the system within the framework of the SBMFT. By noting that the anomalies were also 
observed in experiments for different longitudinal fields at almost the same temperature, 
although with smaller peaks than applying the transverse fields. The fact that our SBMFT 
cannot produce this result shows that there might be some limitations for this approach. 
With increasing the magnetic field, the temperature T c where the field-induced transition 
occurs, moves to low temperature side. Combining the results of susceptibility (see below), 
we may identify the transition as one from a disordered phase to a spin-polarized phase. We 
could understand this phenomenon as a consequence of the broken of the Z 2 symmetry. 16 
When the external field is applied along the z-axis, the XY symmetry is retained. On the 
contrary, when the magnetic field is applied perpendicularly to the z-axis, the XY symmetry 
is broken and an Ising anisotropy is produced (See, e.g. Ref. [8] and references therein). At 
a given field an anomaly of specific heat will thus appear at a given temperature due to the 
Z2 symmetry breaking, which is an indicative of a field-induced transition from a disordered 
phase to a spin-polarized phase as the transverse field tends to align spins perpendicular to 
the easy axis. Here we would like to point out that the anomaly we observed is different 
from the usual broad peak of C v which appears at temperature T p whatever the magnetic 
field is present or absent, while T p decreases with increasing the field. Since the limitation of 
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the approximation we cannot find such a broad peak of C v within the frame of SBMFT, as 
discussed in a number of references. This may be a shortcoming of the SBMFT. However, 
at high temperature we have observed the specific heat decreases with increasing the field, 
which shows a correct trend for S — 1 AFM chain. On the other hand, we note that the 
anomalies in the curves of C v were observed experimentally for NDMAP at high magnetic 
fields. Our calculated result is in qualitatively agreement with this experimental obervation 8 
if we make a rescale for the field as h x — > (h c — h x ), 19 where h c is a critical field at which 
the anomaly of specific heat disappears at zero temperature. In the present case h c ~ 0.90. 
With this field rescaling, we find that the temperature where the anamoly occurs increases 
with the field, consistent with the experimental observation. We would like to point out that 
we could not compare the experimental data directly with our calculated results, because the 
experimental data 8 include the contribution of the lattice which is not given. However, one 
may see that the results obtained on the basis of the SBMFT in the present fashion might 
capture qualitatively some experimental features of the S = 1 HAF chain in the presence of 
a transverse field. Besides, our SBMFT result seems also to be in qualitatively agreement 
with the upper phase boundary of the exact diagonalization of finite chains, although the 
latter results are extracted from the staggered susceptibility. 16 

We have also investigated the effect of the single-ion anisotropy on the specific heat in 
a given field h x = 0.3, as shown in Fig. 2. The result shows that the anisotropy does not 
have so much effect on the specific heat, and it only causes slight changes in the peaks 
and dips of C v . It appears that C v increases, though small, with increasing D, as indicated 
in the inset of Fig. 2. One may notice that this result is very different from the case in 
the presence of a longitudinal applied field, where the effect of the single-ion anisotropy 
on specific heat is various. 10 This result is readily understandable, because the single-ion 
anisotropy is XY-like, and if the magnetic field is applied perpendicular to the z-axis, the 
single-ion anisotropy would have no much effect on the physical obervables, as illustrated 
in the present case, while such an anisotropy would matter if the magnetic field is applied 
along the z-axis, as manifested in Ref. [10]. 

The external field dependence of the specific heat for different temperatures at D — 1.0 
is depicted in Fig. 3. It is shown that at a given temperature the specific heat increases 
with increasing field, and after arriving at a maximum, decreases faster and then slow with 
increasing field. One may see again that there are anomalies in the curves of the specific heat 
versus the external field. It is interesting to note that the magnetic field at which the specific 
heat shows a peak at a given temperature is consistent with those found in Fig. 1. For a given 
temperature, C v (h x ) behaves as C v (h x ) ~ c (T) + ci(T)h l J 2 + c 2 (T)h x + c 3 (T)h%/ 2 , where q 
(i = 0, 1,2,3 ) are temperature-dependent coefficients. For example, at T = 0.5 the fitting 
results give c = 0.1643, a = 0.7013, c 2 = -1.2047, c 3 = 0.468. At h x — ► 0, C(T = 0.5) 
is 0.1643. It can be seen that after making a field recaling h x — > h c — h x one may find 
that the position of the anomaly of the specific heat moves to high field side with increasing 
temperature. 



B. Transverse Susceptibility 

The temperature dependence of susceptibility (x) at D = 1.0 in different fields is depicted 
in Fig. 4. It is observed that the susceptibility shows a broad peak at low temperature in 
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low fields (in our case h x < 0.485), which is a characteristic of S — 1 Heisenberg AFM chain, 
and it goes to zero at T — > for the field less than h c \ = 0.485, suggesting that the system is 
in a Haldane gapped phase in this regime. For 0.485 < h x < 0.9, we found the susceptibility 
of this ID spin system goes to a finite value as T tends to zero, implying that in this regime 
the system is in a Luttinger liquid (LL) phase. It is clear that there exists a transition from 
the Haldane gapped phase to the LL phase at h x = h c \ at which the excited gap closes. 
The result agrees with that of the effective-field theory and size scaling analysis 18 with the 
finite chain calculation 19,20 , which indicates that the field-induced transition occurs from 
the disordered ground state to the gapless Tomonaga-Luttinger liquid phase at some critical 
external field h ci . In addition, a phase boundary in field (/i c i) -temperature (T) plane is thus 
obtained by observing the variation of susceptibility. On the other hand, for h c \ < h < h c2 , 
the gap vanishes and we find a range of linear T dependence for the specific heat, which is also 
a characteristic for the LL phase. 21,22 For the higher field (h x > h c2 = 0.9) the susceptibility 
diverges and persists up to saturation as temperature tends to zero, showing that the system 
now enters into a spin-polarized phase. Thus, there must be a transition from the LL phase 
to a spin-polarized phase at h x = h c2 . However, it is not easy to determine h c2 according 
to the divergence of the susceptibility. Fortunately, we observed that there is a curvature 
change in the curves of x versus temperature for different fields, as enlarged in the inset 
of Fig. 4. Although the curvature change in the susceptibility itself is not an indication 
of a phase transition, by observing that the susceptibility shows different behaviors at low 
temperatures in different phases, one could identify when the system enters into the spin- 
polarized phase from the LL phase by observing the curvature change of the susceptibility. 
The positions of curvature change are the same as the positions in curves of the specific heat 
at which the specific heat exhibits an anomaly, which gives the estimation of h c2 . In other 
words, at h c2 the specific heat shows anomalies, indicating a transition occurs at h c2 . Note 
that the curvature-changing position moves to low temperature side with increasing field. 
It can be observed that the behavior of susceptibility in the transverse field is very different 
from the results in the longitudinal field, 10 where there is no curvature-changing observed 
in the curves of x versus temperature. 

The effect of the single-ion anisotropy on susceptibility is also investigated. It is found 
that the single-ion anisotropy does not change so much the shape and the magnitude of 
susceptibility. Though the effect on susceptibility is not obvious, it can cause a small de- 
creasing of x at a given temperature and a given field. One can obtain a schematic phase 
diagram in field (h x ) -anisotropy (D) plane by observing the variation of susceptibility, as 
shown in Fig. 5. The phase boundary is characterized by closure of the excited gap. The 
system is in a gapless phase above the boundary and is in a gapped Haldane phase below the 
boundary. One may note that the change of h x with D in the diagram is only a few percent. 
Whatever the anisotropy is strong or weak, the system is always in the Haldane phase if 
h x is less than 0.45. When D is larger, the external field required to drive the system into 
the gapless phase becomes larger. For example, at D — 0, the critical field is 0.45, and at 
D = 1.6, it becomes 0.50. It is found that the critical field in the case of the transverse 
field is larger than that in the longitudinal field. For instance, at D = and T — > the 
critical field in the transverse case is 0.45, while it is 0.155 in the longitudinal case. Our 
result is in agreement with the estimation on the critical field at D = on the basis of 
exact diagonalization of finite chains, 16 and is qualitatively consistent with the experimental 
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extrapolation in NDMAP 8 . 

Fig. 6 gives the field-dependence of susceptibility at D — 1.0 at different temperatures. 
It is observed that for h x < 0.9 the behavior of x versus transverse field at low temperature 
is quite different from that at high temperature. The positions of curvature change in x are 
consistent with our observation for the anomalies in specific heat. We also observed that the 
single-ion anisotropy does not have much effect on the behavior of the susceptibility versus 
transverse field. 

III. SUMMARY AND DISCUSSION 

By summarizing the above results on the specific heat and the susceptibility, we may 
propose a schematic phase diagram in field (/i^-temperature (T) plane within the framework 
of SBMFT. Since the anisotropy does not have much effect on the behavior of thermodynamic 
obervables of the system under interest, without loss of generality we present the phase 
diagram for D = as an example, as shown in Fig. 7. When the transverse field is less 
than h c i, the system should be in a gapped phase (Haldane phase). At h x = h c \ the gap 
closes. The lower boundary for h c i is determined by observing the closure of spin gap in the 
curves of the susceptibility as a function of the applied transverse field and temperature. For 
< h x < h C 2, the system might be in a Luttinger liquid phase, characterized by finite values of 
the susceptibility at T — > 0. When h x > h C 2, the system goes into a spin-polarized phase. 
At h x = h C 2, the anomalies in the curves of the specific heat as well as the susceptibility 
appear. It is interesting to observe that the shape of the LL phase becomes more symmetric 
for D nonzero. It seems that our proposed phase diagram based on the SBMFT is somewhat 
different from the results of exact diagonalization on finite chains. 16 In the phase diagram 
presented in Ref. [16] the lower and upper boundaries of the antiferromagnetic ordered 
phase are given by h c \ and h C 2, which are determined by the staggered susceptibility. While 
in our proposed phase diagram, h c \ and h c2 are determined by observing the behaviors 
of the uniform transverse susceptibility as well as the specific heat. Considering the mean- 
field feature of our method, we cannot determine with accuracy whether a small region of an 
ordered phase at very low temperatures (T < 0.1) in the LL phase exists, as the susceptibility 
we calculated is not a staggered but uniform one. In this sense, our proposed phase diagram 
is not incompartible with that of Ref. [16], which can be viewed as a complement for finite- 
size calculations. Therefore, it could be better to understand the thermodynamic properties 
of this ID spin system in the presence of the applied field based on a combination of our 
SBMFT results and finite-size calculations. 
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APPENDIX A: FORMALISM OF THE SBMFT 



In this Appendix, we shall present the detail derivation of the self-consistent equations 
based on the SBMFT. These equations are applied to get the thermodynamic properties of 
the system. 

We suppose that the system in Eq.(l) with J > is defined on a bipartite lattice with 
sublattices A and B. For sublattice A, we denote the spin operator Sf by Schwinger bosons 
cij and bi (see Ref. [11]): 

<7 A — —n + h + <7 A — —nh 

5^ = ^ + 6^ + 1), teA {A1) 

satisfying the constraint Sf = (afdi + bfbi) for each site on sublattice A, and for sublattice 
B 

qB _ _ L. CB _ n + U+ 

with = (alfaj + Vjbj) for each site on sublattice B. The bosons {a} and {b} obey the 
standard commutation relations. Although one may note that the definitions in Eqs. (2) and 
(3) introduced in Ref. [11] are slightly different from the standard form, we find that such 
a form is quite convenient for our purpose. On account of the definitions with constraints, 
the Hamiltonian (1) for S — 1 can be rewritten as 

H = H + H D + H ex , (A3) 

(A4) 



Ho — JJ2(i,j) Sj • Sj — —2JJ2i=l,8 ^-i,i+S^i,i+S 
~JJ2i=l,6 a i 0'iO'i+S a i+S + J J2i=l,8 a i a ii 



3D 

Hd = D £(S*) 2 = — ]T(a+a, + bfh) + DN/2, (A5) 
i=i 4 i=i 

H ex = -lib Y\9dSih x + Sfh y ) + 9ll S*h z ], (A6) 

i 

where A iji+ s = ^(afbf +s + bidi + s), and 2N is the total number of sites. To implement the 
constraints J2ieA( a t a i + K^d = S A an d J2jeB( a ^ a j + — $ B where S A ^ denotes the 

total spin of sublattice A(B), we should introduce two kinds of Lagrangian multipliers \f 
and into the system. At the mean-field level, for the sake of simplicity we may take the 
average value of the bond operator < A i>i+ s >= A to be uniform and static, and so are 
< \f >= \ A and < Af >= \ B . 

The mean-field Hamiltonian in momentum space reads 

H = H + H D + H ex , (A7) 
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Ho = -JE k {A*(at k bt + b k a_ k )zr k + (b k a- k + a + _ k b+)z lk A} 
-JJ2k(< n i+ s,a > at k a- k + < n iA > at k a^ k ) 

+ JE fc a± k a- k + \ A £ fc [(at fc a_ fc + b + k b k ) - S A ] (A8) 
+\ B E k [(at k a- k + bth) -S B ] 
+2JzNA*A + JNz < n ija >< n i+ha >, 



ID 

Hd = — E( a + Kh) + DN/2, (A9) 



1 — 1 

Hex = T^f^B ^2[-g±h x (at k bt + b k a- k ) + —g±.h y {at k bt - a- k b k ) + g\\h z (at k a^ k + b\b k + 1)], 

1 k 3 

(A10) 

where z is the number of nearest neighbor sites, j is a complex number, and r ) k — \ J2s elkS — 
cos k. The sum over fc is restricted to the reduced first Brillouin zone. 
Utilizing the Bogoliubov transformation 

a_fc = cosh9 k a k + smh9 k /3 k , b k = smli^a^ + cosh9 k (3 k , (All) 

with 6 given by 

tanh29 = ~^h x +Je(2JzA lk ) ^ ^ 

with 

2A' = (\ A + \ B ) - B/2 + 3D /A, B = J[(< n m , a > + < n l>a >) - 1)], 

we obtain the energy spectrum of the system in the presence of a transverse external field 
hex = (h x , 0,0): 

E k = v/(2A') 2 - (2it>e \JzA lk \ - fj, B g ± h x y - B/2, (A13) 

El = y/(2A'¥ - (2Re \JzA lk \ - fj, B g ± h x y + B/2. (A14) 

The free energy per site takes the form of 

& = ^£ fe {m[2sinh §(££)] +ln[2sinhf (££)]} 
+JzA*A + B/4 - (S A /2 + 1/2)A A - X B (S B /2 + 1/2) (A15) 
-D/2 + Jz< n ija >< n i+ha > /2. 

The mean-field self-consistent equations can be obtained by minimizing the free energy. 
Without loss of generality, we may set X A = X B . Then, we perform 5F/5X A = 0, 5F/5A* = 
0, 8F/8 < n i+s ,a >= and 5F/5 < n iA >= 0. Rescale the parameters (A A , \ B ,A, (3) 
— ► (A, V, «) : 
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JAz n Ak . . , „. 

77 = — - , (3 = — . (A16) 

A z 

Then, the angle in the Bogoliubov transformation can be expressed in a compact form 

1 _fJ-B9±h x , 

cosh2fl fc = -p sinh2fl fc = -= 2A ' (A17) 
The self-consistent equations become 



+ S*)/2 + 1 = ^ E fc -7 i— ={coth «(2AVl - (777, - ^) 2 + 5/2) 



+ coth K (2Ayi - ( Vlk - ™^) 2 - 5/2)}, 



(A18) 



A' = £ E,*> ■ -7^^7 ^° th «(2AVl-(^"^) a + 5/2) 



+ coth K (2A V / 1 - (777. - ^") 2 - 5/2)}, 



4 < n i>0 > +1 = 

^E*{H - 1] coth«(2Ayi - (777, - ™^) 2 + 5/2) 

Vi-Cto-^^t^) 2 (A20) 

+[- 1— - - 1] coth «(2A7l - (777, - ^) 2 - 5/2)}, 

4 < Wi+La > +1 = 

* E * V-^-^ + 1] coth ^Vi-(^-W - */') (A21) 

+[- L_= - 1] coth K (2A v /l-(777,-^) 2 + 5/2)}. 

As is easily seen, the aforementioned equations are highly coupled, and it is not possible to 
get useful analytic results. They are only solved numerically to get the relavant thermody- 
namic quantities. 
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FIGURE CAPTIONS 



Fig. 1 The temperature dependence of the specific heat at D = for different transverse 
fields. 

Fig. 2 The temperature dependence of the specific heat for different anisotropies at 
h x = 0.3. Inset: The enlarged part of the specific heat versus temperature for different 
anisotropies at h x = 0.3. 

Fig. 3 The ft^-dependence of the specific heat at D — 1.0 for different temperatures. 

Fig. 4 The temperature dependence of the susceptibility at D — 1.0 for different trans- 
verse fields. Inset: The enlarged part of the susceptibility versus temperature at D = 1.0 
for different transverse fields. 

Fig. 5 Schematic phase diagram in h x — D plane. The solid line is the phase boundary 
on which the gap closes. 

Fig. 6 The /^-dependence of the susceptibility at D = 1.0 for different temperatures. 

Fig. 7 Schematic phase diagram in h x — T plane. 
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